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A detailed theory describing linear optics of vapors comprised of interacting multi-level quantum 
emitters is proposed. It is shown both by direct integration of Maxwell-Bloch equations and using 
a simple analytical model that at large densities narrow transparency windows appear in otherwise 
completely opaque spectra. The existence of such windows is attributed to overlapping resonances. 
This effect, first introduced for three-level systems in [R. Puthumpally-Joseph, M. Sukharev, O. 
Atabek and E. Charron, Phys. Rev. Lett. 113 , 163603 (2014)], is due to strongly enhanced dipole- 
dipole interactions at high emitters’ densities. The presented theory extends this effect to the case 
of multilevel systems. The theory is applied to the D1 transitions of interacting ®®Rb atoms. It 
is shown that at high atomic densities, ®®Rb atoms can behave as three-level emitters exhibiting 
all the properties of dipole induced electromagnetic transparency. Applications including slow light 
and laser pulse shaping are also proposed. 
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I. INTRODUCTION 

The research of phenomena related to light-matter in¬ 
teraction has been experiencing signihcant growth for the 
past decade [1] in the domain of nano-optics. From the 
fundamental point of view, one of the most intriguing 
question lies in the dynamics and in the optical prop¬ 
erties of many-body systems. For instance, new col¬ 
lective modes were predicted and measured in materi¬ 
als composed of interacting quantum emitters (atoms, 
molecules, etc.) [2]. Experiments in that regard are 
noticeably advanced compared to theory, as it has long 
been realized that treating multi-body systems is quite 
challenging. Even systems comprised of just a few inter¬ 
acting quantum dipoles are difficult to investigate [3, 4]. 
Semi-classical descriptions of light-matter interaction, in 
which a classical electric field interacts with a quantum 
system, enormously simplify modeling leading to results 
that both support experiments and predict new phenom¬ 
ena [2, 5, 6]. 

Collective effects [7] arise from the fact that the elec¬ 
tric field experienced by an arbitrary quantum emitter 
is the sum of both the applied held and the radiation 
from all other atoms in the system under consideration. 
When it is strong enough, the radiation from the dipoles 
set them coupled. This is the source of cooperativity. 
One of the important fundamental cooperative effects, 
the collective Lamb shift which is due to the exchange of 
virtual photons between the particles in a dense sample of 
atoms was recently observed [8] . Dependent of the coop¬ 
erativeness in homogeneous and inhomogeneous systems 
were studied both experimentally [9, 10] and theoretically 
[4, 11, 12]. At present, research on those cooperative ef¬ 
fects concerns on both cold [13-17] and thermal [8, 18, 19] 
samples of atoms since they offer the observation of dif¬ 
ferent interesting quantum phenomena. Because of the 
ability of light to carry effectively quantum information. 


understanding the optical response and the light scatter¬ 
ing in dense samples is also interesting for many appli¬ 
cations such as slow and stopped light [20, 21], quantum 
information [22, 23], trapped light [24] and optical mem¬ 
ories [25]. 

In the linear regime, the propagation of light through 
any medium can be investigated by extracting its elec¬ 
tric susceptibility [26, 27]. It was shown recently [6] 
that the electric susceptibility of systems comprised of 
three level quantum emitters at high densities exhibit a 
Fano-type profile that results in a new phenomenon of 
induced transparency called Dipole Induced Electromag¬ 
netic Transparency (DIET). DIET is due to the collective 
response of the system towards the external electromag¬ 
netic (EM) held via dipole-dipole interactions. Following 
this work, in this paper, we present a detailed theory of 
DIET and extend the model to multilevel systems. We 
also apply our theory to a realistic system, namely the 
D1 line of ®5Rb. 

The paper is organized as follows. The next section 
discusses in detail our theoretical description for Al-level 
systems and the analytic model used to analyze the re¬ 
sults. Section III provides results and a comprehensive 
discussion of how linear optical properties of a system 
comprised of interacting two-level emitters depends on 
various material parameters such as the particle den¬ 
sity. In Section IV, we extend our discussion to sys¬ 
tems having more than one excited states. First we con¬ 
sider interacting three-level emitters. The results pro¬ 
vided demonstrate dipole induced electromagnetic trans¬ 
parency (DIET). This phenomenon is explained using an 
analytical model and is extended to the multilevel D1 
transitions of ®^Rb. Section V discusses possible appli¬ 
cations of our results and hnally, we conclude our work 
in Section VI. 
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II. THEORETICAL MODELS 

Our goal is to describe the linear dynamics of quantum 
emitters in a self-consistent manner taking into account 
collective effects. The most comprehensive and thus com¬ 
plete description would be to implement a fully quantum 
model based on the Jaynes-Cummings Hamiltonian [28]. 
On the other hand, the quasi-classical theory of light- 
matter interaction recently applied to hybrid materials 
has proven to lead to both qualitative and quantitative 
agreement with experiments [29]. Keeping also in mind 
that the majority of experiments in nano-optics, such 
as white light spectroscopy, for instance, are performed 
under conditions corresponding to high occupation num¬ 
bers of photons in a given EM mode, the quasi-classical 
description is well justified. 

The idea is to separate the description of EM waves 
from the dynamics of quantum emitters considering the 
spatio-temporal evolution of EM radiation using classical 
Maxwell’s equations while applying the full machinery of 
quantum mechanics to describe the response of quantum 
emitters to EM excitation. This results in a system of 
coupled Maxwell-Bloch equations [30]. If one wants to 
observe collective effects, it is imperative to solve the cor¬ 
responding equations of motion self-consistently without 
any decoupling [31]. 


ground state |0) via the time dependent EM interaction 
Vint {f, t). The density matrix p (r, t) satisfies the follow¬ 
ing Liouville-von Neumann equation 

= (4) 

where H = Ho + Vnt (r, t) is the total Hamiltonian and T 
is the superoperator describing relaxation and dephasing 
processes taken in the Lindblad form under Markov ap¬ 
proximation [32] . Off-diagonal elements of T include the 
pure dephasing rate 7 * and the diagonal elements include 
the radiationless decay rate T of the excited states. 

The field free Hamiltonian Ho is given by 

N-l 

Ho='^f>^j\j){j\- (5) 

3=0 

We define the Bohr frequencies as luoj = uJj — uio and 
the interaction of a single emitter with the EM fields is 
written as 

N-l 

Vint {r,t) = ^^j(^’^)(|0)(j| + b)(0|): (6) 


A. Theoretical Model and Numerical Simulations 


The dynamics of the EM field components, E and Lf, 
is simulated using Maxwell’s equations 


9H „ - 

Po^ = -V X E 
at 

dt dt 


(la) 

(lb) 


where Eq and po ^re the permittivity and the permeabil¬ 
ity of free space, respectively, and P{r, t) is the macro¬ 
scopic polarization of the system at position r and time 
t. The dynamics of the latter is considered both quan¬ 
tum mechanically and classically. In both cases we as¬ 
sume that the emitters are continuously distributed in 
space and we neglect static emitter-emitter interactions. 
Under such conditions one can express the macroscopic 
polarization as 


P{f,t) = no{p), 


( 2 ) 


where no is the number density of emitters. The average 
dipole moment, {p), is either determined directly from 
classical equations of motion or quantum mechanically 
by evaluating 

{p) =TT:[p{r,t)p], (3) 

where p(r, t) denotes the density matrix of the system. 
Each quantum emitter has {N — 1) excited states rep¬ 
resented by \j), with j > 1 and N > 2 , coupled to the 


where {Hj(r, t)} are the instantaneous Rabi frequencies 
associated with the interaction between the quantum sys¬ 
tem and the local EM fields. In the present work we as¬ 
sume that the excited states are not directly coupled to 
each other but only to the ground state |0). The corre¬ 
sponding system of coupled differential equations reads 
[33] 


Poo 

P33 

P 03 


- P30) + ^Pjj, 


J>1 

t){pjo - Poj) - Tpjj, 


zUj(r,t)(poo 


P33) + 


iuoj - 7 


POj) 


(7a) 

(7b) 

(7c) 


where 


7 = 


27 * -bT 
2 


( 8 ) 


Eqs. (1) and (4), coupled through the evolution of the 
quantum polarization ( 2 ), form the basis of the model. 
As this system of equations is propagated in space and 
time on a grid, one should note that each grid point is 
effectively a point-wise dipole with the amplitude of an 
individual emitter times the number density. This es¬ 
sentially means that dipole-dipole interactions within a 
single grid point are neglected. This approximation al¬ 
though valid at low densities has to be corrected for high 
densities. Treating the dipoles exactly to include those 
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contributions is extremely difficult and it is almost im¬ 
possible to solve the system of equations by using cur¬ 
rent computational facilities [4]. An alternative way to 
include dipole-dipole interactions at a single grid point 
level is to introduce the well-known mean field Lorentz- 
Lorenz correction term to the local field [26], with 

Aiocal = E + -—. (9) 

3^0 

This local electric held enters the Liouville-von 
Neumann equation (7) through the Rabi frequencies 
= iJLojE\oca\(f,t)/fi, where /roj denotes the transi¬ 
tion dipole moment between states | 0 ) and |j). 

Using the proposed model we simulate the linear opti¬ 
cal response of a vapor comprised of interacting quantum 
emitters as schematically depicted in Fig. 1. To simulate 
a hrst order elastic scattering/absorption we implement 
the short pulse method [30]. A weak ultra-short incident 
pulse polarized along x is launched at normal incidence 
{z direction) on a slab of quantum emitters of hnite thick¬ 
ness. 



FIG. 1. (Color online) The schematic setup of simulations: a 
slab of thickness £ composed of interacting quantum emitters 
is exposed to linearly polarized EM fields acting normal to 
the interface. 


For such a geometry, Maxwell’s equations (1) reduce 
to 


dHy _ dE, 
dt dz ’ 

^ dE^ _ dHy dP^ 
~ ~lfz W' 


( 10 a) 

( 10 b) 


Maxwell’s equations are numerically integrated using 
a hnite-difference time-domain (FDTD) technique [34] 
with a spatial grid step, Sz = 1 nm and temporal grid 
step, St = 1.67 as. Convolutional perfectly matched 
layers (CPML) [35] absorbing boundaries of thickness 
19 nm are used to avoid non-physical reflections from 
the boundaries of the simulation region. Concurrently, 
Liouville-von Newmann’s equations (7) are propagated 


in time using the fourth order Runge-Kutta method, as¬ 
suming that all emitters are initially in their ground state 
|0). Light excitation is provided by a gaussian pulse of 
FWHM Suj = 27r/T, where r is the pulse duration. 

To evaluate transmission, reflection and extinction, we 
calculate the EM energy flux S{uj) on the input and the 
output sides of the slab using 


Siuj) 


ExHy 


\Ex^incEy^ii 


( 11 ) 


where E^, Hy, Ex^mc and i7y,inc are the Fourier compo¬ 
nents of the total and incident EM helds. Eq. (11) gives 
the reflection spectrum R{uj) if Ex and Hy are the re¬ 
flected EM fields and the transmission spectrum T(w) if 
Ex and Hy are the transmitted EM helds. Both spectra 
carry the information about the dynamics of the system. 

And hnally the extinction A{uj) is estimated as 


A{uj) = 1 - [T{uj) + R{uj)]. (12) 

The layer of quantum emitters is placed at a distance 
of 1.25 ^m from the source. The reflected held is mea¬ 
sured at a point between the source and the layer which 
is at 10 nm away from the source and the transmitted 
held is measured on the output side, 2.5 /rm far from the 
source. We verihed that the results are independent of 
those distances as long as the pulse acts normal to the 
layer. 


B. Semi-classical Approximation 

In the linear regime, poo lor J ^ ^od the 

homogeneous system of A^-level quantum emitters can 
be described by a set of {N —1) coupled harmonic dipole 
oscillators. The macroscopic polarization associated with 
the transition between | 0 ) and |j), Pj, follows a classical 
equation of motion 


( 13 ) 

where 7 cz = 27 and ruj is the effective mass of the os¬ 
cillating charge Qj driven by the electric held i?iocai(i)- 
Taking into account Eq. (9) and subsequently applying 
a Fourier transform, we obtain 


W,{uj)P,{iz) = eoiz^.Exiuj) + ^ Pk{uj), (14) 


where 


Wj{u}) = uj^y - Uj'^ + ijciuj (15) 
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and where uj-pj is the plasma frequency representing the 
collective oscillations of the dipoles in the system 


LU 


2 _ 
■Pj - 


noq^j 

eorrij' 


(16) 


With the assumption that the maximum amplitude of 
dipole oscillations in the absence of a driving field is given 
by the quantum harmonic oscillation length, we have 


2 _ 2a;ojnoAij 

Sfeo 


(17) 


In weak electric fields where Pj{u>) = XjSoExi^^), Eq. 
(14) gives the susceptibility associated with the 

j*** dipole. Summing over all dipole oscillators in the 
material, we obtain the total susceptibility of the system 
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(18a) 


III. INTERACTING TWO-LEVEL EMITTERS 

First we consider a slab of interacting simple two-level 
emitters characterized by the transition energy fno^i = 2 
eV, the transition dipole = 2 D, and the total deco¬ 
herence rate 7 = 10.5 THz. We calculate the response of 
the system to an extremely short pulse of width t = 0.18 
fs using coupled Maxwell-Liouville-von Neumann equa¬ 
tions in different collective interaction regimes [36]. The 
results are analyzed using the semi-classical approxima¬ 
tion described in Section IIB. 


A. Semi-classical approximation and the Spectra 

The electric susceptibility of interacting two-level emit¬ 
ters in the semi-classical model can be written as 


X{^) 


, ,2 


(19) 


where woi = is the shifted resonant fre¬ 

quency of the system. The shift in transition frequency, 
Ai, which is due to the Lorentz-Lorenz correction in¬ 
cluded in the formulation of the local electric field via 
Eq. (9) can be estimated from the semi-classical model. 
The first order term gives the well-known Lorentz-Lorenz 
shift (LL shift) [12, 26, 37]. 


( 20 ) 


The case of strong dipole-dipole interactions is char¬ 
acterized by Ai 7 , while Ai <C 7 corresponds to the 
weak interaction regime. 

Typical reflection spectra are shown in Fig. 2 as a 
function of the reduced detuning, d = (w — a;oi)/ 7 , for 
different values of A 1/7 for a slab of thickness t = 600 
nm. In the weak interaction regime (A 1/7 = 10“^), the 
system reflects a very small amount of the total energy (of 
the order of 10“®) around the transition frequency. Such 
a response is similar to that of a single emitter since it 
is characterized by a typical Lorentzian profile. As the 
value of A 1/7 increases, and especially for Ai > 7 , the 
interaction between the dipoles becomes dominant and 
the system starts to respond to the applied field collec¬ 
tively. For A 1/7 > 10, the reflection spectrum broad¬ 
ens resulting in a window of frequencies within which 
more than 80% of incident energy is reflected. This re¬ 
flection window was explained in 2000 by R. J. Glauber 
and S. Prasad in Ref. [11, 12] by considering the ex¬ 
ponentially decaying modes of coherent excitation of the 
medium that depend strongly on the frequency. 

For a slab of thickness I larger than the dipole wave¬ 
length Aoi = 27rc/woi) the reflection at the vacuum/slab 
interface (shown as interface 1 in Fig. 1) dominates and 
the contribution from the possible multiple reflections in¬ 
side the slab can be neglected. For such a slab, the width 
of the reflection window can be estimated from the reflec¬ 
tion probability at this interface. 


7^(a;) 


1 - M{uj) " 
1+A7(a;) ’ 


( 21 ) 


where A/’(a;) = Re[y^l -I- A(a;)] is the real part of the 
effective refractive index of the system. In the case 
of total reflection, 7^(w) = 1 and the real part of the 
electric susceptibility is —1. Equating the real part of 
Eq. (19) to —1 leads to the range of frequencies within 
which the reflection probability reaches its maximum: 
Woi — Ai < w < Wqi -b 2 Ai and hence the frequency 
window with the maximum reflection shows a width of 
3Ai which agrees with the quantum simulation [ 6 ]. For 
relatively thin slab of width I < Aqi, the width of the 
reflection window is noticeably below 3Ai. This is due 
to the fact that in this analysis we have neglected con¬ 
tributions from multiple reflections that occur inside the 
slab. For a thin system, contributions from the two dif¬ 
ferent vacuum-medium interfaces are not negligible and 
obviously affect the reflection spectra. 

Fig. 3 shows the extinction (panel (a)) and transmis¬ 
sion (panel (b)) spectra at A 1/7 = 10“® (blue solid line) 
and A 1/7 = 12 (green dashed line). For A 1/7 = 10“®, 
almost all the incident energy is transmitted. Clearly in 
the weak interaction regime, the local polarization barely 
modifies the local electric field experienced by a given 
emitter. Hence the emitters respond to the field as if 
they were independent. As the mutual interaction be¬ 
tween the dipoles increases, the system becomes opaque 


Ai = noM?/9feo, 
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FIG. 2. (Color online) Reflection spectra as a function of 
the reduced detuning S = (ui — uJoi)/j for a slab of interact¬ 
ing two-level emitters in different interaction regimes. The 
thickness is ^ = 600 nm and the transition dipole /ii = 2 D. 
The radiationless decay rate and the pure dephasing rate are 
r = 0.5 THz and 7 * = 10 THz , respectively. Reflection for 
small values of A 1/7 is weak and hence amplified (for better 
comparison) by the constant factors shown in the Figure. 

and reflects most of the incident energy [12]. In com¬ 
plement with reflection and transmission, the remaining 
energy is absorbed by the system and decays due to ra¬ 
diationless transitions described by F in Eq. (4) . In par¬ 
ticular at A 1/7 = 12, the system loses its transparency, 
as seen in Fig. 3(b), and behaves almost like a mirror 
around a broad range of frequencies near the transition 
frequency. 



FIG. 3. (Color online) Extinction (panel a) and transmission 
(panel (b)) probabilities as functions of the reduced detun¬ 
ing S. The solid blue curves (left vertical scale) represent the 
results for A 1/7 = 10~^ and the dotted green lines (right 
vertical scale) show the data for A 1/7 = 12. All other pa¬ 
rameters are the same as in the Figure 2. 


It should be emphasized that the reflection from 
strongly interacting dipoles resembles Bragg mirrors (also 
referred to in the literature as distributed Bragg reflec¬ 
tors) [38] . Bragg mirrors are composed of alternating di¬ 
electric layers with different refractive indices. The high 


reflectivity is achieved by multiple reflection from layers 
of the mirror. This leads to a controllable constructive 
interference. Spatially varied refractive index profiles de¬ 
termine the fringes in the spectrum. In contrast to Bragg 
mirrors our system is completely homogeneous. Since we 
perform simulations in a linear regime, the refractive in¬ 
dex of the system is a constant throughout the layer. 
This rules out the possibility of having multilayer inter¬ 
ferences. All the effects observed in the spectra are due 
to the collective response of the system. The width and 
the position of the reflection window are characterized by 
the physical properties of the system such as the transi¬ 
tion dipole, transition frequency, material density. The 
qualitative explanation for the modihcations in the spec¬ 
tra when the system enters the strong interaction regime 
relies on the fact that the polarization of the system is 
enhanced due to the strong dipole-dipole interactions. 
Prior to the detailed quantitative analysis of the spectral 
features, we want to estimate the accuracy of the semi- 
classical model and compare it with results obtained via 
exact numerical simulations. 

Fig. 4 shows the relative errors of the electric suscep¬ 
tibilities for interacting two-level emitters obtained using 
the semi-classical model with respect to the numerical in¬ 
tegration of the Maxwell-Bloch equations, Eq. (19) in the 
weak (A 1/7 = 10“^ ) and strong (A 1/7 = 12) interac¬ 
tion regimes. The solid blue line shows the relative error 
for weakly interacting dipoles and the red dashed line, is 
the same for strongly interacting dipoles. Both relative 
errors are peaked around the corresponding transition 
frequencies with a maximum amplitude of about 0.08% 
and we conclude that the analytical model is accurate 
enough to analyze the numerical results. 



FIG. 4. (Golor online) Relative error in the absolute val¬ 
ues of the electric susceptibility using the semi-classical ap¬ 
proximation with respect to the numerical integration of the 
Maxwell-Bloch equations as a function of the reduced detun¬ 
ing for parameters used in Figure 3. The blue solid line shows 
the relative error for weakly interacting dipoles and the red 
dashed line presents the same for strongly interacting dipoles. 
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B. Fabry-Perot Modes 

The semi-classical model explains the spectral features 
such as the shift in the resonance frequency and the width 
of the reflection window for a system of interacting two- 
level emitters in the strong interaction regime. The fea¬ 
tures that are not yet interpreted by the model are side¬ 
bands (see Fig. 2 for A 1/7 = 10). Those however can 
be interpreted as simple Fabry-Perot modes [ 8 , 39]. The 
semi-classical model does not yet contain any informa¬ 
tion about the geometry of our system, assuming homo¬ 
geneous distribution of the emitters. With two or more 
interfaces present, one can expect to observe interferences 
between reflected and transmitted fields due to the pos¬ 
sible multiple reflections inside the slab. In this case, the 
system can be seen as a Fabry-Perot etalon consisting of 
two parallel partially reflecting interfaces shown as the 
interfaces 1 and 2 in Fig. 1. EM helds at frequencies 
outside the transmission window are partially reflected 
and transmitted, resulting in non-zero optical path dif¬ 
ferences and hence giving rise to an interference pattern 
as seen in Fig. 5 for instance. Introducing the decay of 
the EM helds inside the slab and taking into account mul¬ 
tiple rehections with proper boundary conditions [39] , we 
hnd the rehection r{uj) and transmission t(ui) coefficients 
for the system as 

— 2 isin(n(a;)fc£)[n(a;)^ — 

^ [n{uj) + 1]2 - [n(a;) - l] 2 e-* 2 «Mfc^’ ^ 

4n(uj)e~^^ /or,i N 

^ [n{uj) + 1]2 - [n(a;) - l] 2 e-* 2 n(^)fc^ > ^ > 

where n(u}) = \/l + X{uj) is the refractive index, k = 
2tt/X is the propagation constant of the applied held in 
the vacuum and n is the imaginary part of the propa¬ 
gation constant inside the medium. The associated re¬ 
hection and transmission spectra predicted by the semi- 
classical model are given by 

i?(cc) = (23a) 

r(w) = \t{uj)\^. (23b) 

The expressions (23) together with (22) can accurately 
reproduce the positions of the maxima and minima of 
the sidebands in the calculated spectra. Fig. 5 shows the 
rehection spectrum calculated by integrating Maxwell- 
Bloch equations (blue dashed line) and the same obtained 
from the equation (23 a) shown by the solid red line. The 
two curves are in perfect agreement. 

IV. MULTILEVEL SYSTEMS 

Let us consider a ground state coupled to (iV — 1) ex¬ 
cited states. As in the previous case, the inhnite slab 
of emitters is exposed to a transverse electric held. In 



FIG. 5. (Color online) Reflection as a function of the de¬ 
tuning for £ = 600 nm . The dashed blue line is the reflec¬ 
tion probability calculated by integrating the Maxwell-Bloch 
equations and the solid red line is the same obtained from the 
semi-classical approximation (Eq. (23 a)). Other parameters 
are A 1/7 = 12, /ii = 2 D and 7 = 10.5 THz . 


the strong interaction regime, the Lorentz-Lorenz cor¬ 
rection results in redshifts for all transitions. We thus 
expect to observe {N — 1) rehection windows of width 
3Aj each. The most interesting case is when the transi¬ 
tions are signihcantly overlapping, leading to interference 
effects. For three level systems, in the presence of such 
overlapping resonances, we observed DIET which is due 
to the interference between two indistinguishable excita¬ 
tion pathways [ 6 ]. 

The calculated rehection spectrum for a three level 
system with overlapping transitions is shown in Fig. 6 
as a blue solid line. It is compared with the rehection 
spectrum (red dotted line) obtained from a slab of two 
types of uncoupled two-level emitters having the same 
physical parameters which is calculated as the product 
of the rehection probabilities of each dipole. For the 
uncoupled dipoles the reflection spectrum behaves com¬ 
pletely differently compared to the case when they are 
coupled. For uncoupled dipole systems, the total reflec¬ 
tion is peaked where the two independent reflection win¬ 
dows corresponding to each dipole are contributing. This 
region corresponds to a clear overlap of two resonances. 

The strong interaction between the dipoles in the sys¬ 
tem broadens the transition that results in the overlap¬ 
ping of two closely spaced allowed transitions and hence 
sets up a competition between two possible excitation 
pathways. We obtain a narrow window where the ex¬ 
citation probability becomes zero. The presence of a 
hole in otherwise flat reflection spectra is a clear signa¬ 
ture of the destructive interference between the radia¬ 
tion emitted by coupled dipoles and hence a signature of 
DIET [ 6 ]. Except for the observed hole in the reflection 
spectrum, all features discussed in the previous section 
are clearly present in the three-level system spectra. At 
the frequency where there is no reflection, we observe a 
strong transmission if damping processes in the system 
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FIG. 6 . (Color online) The blue solid line is the reflection 
as a function of the reduced detuning of three-level emitters 
with A 1/7 = A 2/7 = 12. The red dotted line is the reflection 
from the same system in the absence of coupling between 
the dipoles. The excited states are separated by a ;2 — oJi = 
4 . 67 . The reduced detuning 5, is defined with respect to the 
transition frequency of the first excited state of this three level 
system. 


are not too efficient. The observed transparency, sim¬ 
ilar to electromagnetically induced transparency (EIT) 
[40, 41], may lead to many potential applications such as 
slow light discussed later in this paper. The position of 
the transparency window can be controlled to a certain 
extend by changing the material parameters [ 6 ]. 

The same type of interference effect can happen if the 
system has a series of closely spaced energy levels. One 
of the potential actual systems for studying DIET are 
Rb atoms. Rb is used experimentally both as a thermal 
gas [ 8 ] and as a cold atomic gas [10]. Dipole-dipole in¬ 
teraction depends on the square of the transition dipoles 
and Rb atoms are characterized by large S to P tran¬ 
sition dipoles [42]. We consider the D1 transition, i.e, 
5 ^ 5 'i /2 to 5 ^Pi /2 of ®®Rb. Both 5 ^ 5 'i /2 and 5 ^Pi /2 are 
split into two sub-levels due to the hyperfine interaction 
[43] and hence there are four dipole-allowed transitions. 
Thus the D1 transitions of ®®Rb isotope can be consid¬ 
ered as a superposition of two three-level systems having 
same excited states but two different ground states. Since 
the energy splitting between the excited states is smaller 
than the splitting of the ground states, it gives a possi¬ 
bility to observe DIET between two excited states that 
are coupled to a specific ground state as well as between 
two excited states coupled to two different ground states 
for different atomic densities. 

The susceptibility of a system consisting of two ground 
states can be written as 


written as 


Aff' — AqSff'j (25) 

where Aq = no/iQ/9feo is the LL shift defined in the 
absence of hyperfine splitting and Sff' is a measure of 
relative strength of the transitions [42]. 

The reflection and susceptibility of a slab of thickness 
£ = 600 nm comprised of ®®Rb for A 0/7 = 3.8 x 10“^ are 
shown in Fig. 7. The reduced detuning S is defined with 
respect to 377.107 THz, the D1 transition frequency in 
the absence of hyperfine splitting [42] . 



Reduced Detuning 5 


FIG. 7. (Color online) Susceptibility (panel (a), log scale) 
and reflection probability (panel (b), linear scale) as a func¬ 
tion of the reduced detuning with A 0/7 = 3.8 x 10“® . The 
reduced detuning 5 is defined with respect to the transition 
frequency of 5 ^ 5 i /2 to 5^Pi/2 in the absence of hyperfine split¬ 
ting (377.107 THz) and 7 = 36.13 MHz. The positions of 
transitions F — 2,3 ^ F' = 2,3 are indicated by the red 
dashed lines marked with {F, F'). 


A(w) = 


+ Thy /^o' 


i-KE 


■3 


2 jW, 


- Ty ^V3‘ 


TWyY 


(24) 


where the index j corresponds to the transitions from the 
sub-level F = 2 and j' is for the transitions from F = 3. 
It is important to note that Eq. (24) can also be used to 
analyze systems composed of different types of atoms or 
molecules. The LL shifts for D1 transitions of Rb can be 


The double peaks in the negative detuning region are 
due to the transitions F = 3 —)■ F' = 2 and F = 3 —>■ 
F' = 3 and those in the positive detuning region are due 
to the transitions F = 2 —F' = 2 and F = 2 ^ F' — 3. 
The separation of the double peaks is due to the hyperfine 
splitting of the ground state. Since the system is in the 
weak interaction regime, the transitions are independent 
and they are not overlapping. It explains the reflection 

























spectrum in Fig. 7(b) with a series of Lorentzian profiles 
corresponding to different transitions. 

Since the dipoles are relatively large for D1 transi¬ 
tions, relatively small increments in the gas pressure will 
couple the different transitions and hence lead to the 
DIET regime. The calculated reflection and transmis¬ 
sion spectra and the susceptibility for such a system with 
^ 0/7 = 21 are given in Fig. 8. Panel (a) shows the 
susceptibility, panel (b) presents the reflection spectrum, 
and panel (c) shows the transmission. 



Reduced Detuning S 


FIG. 8 . (Color online) Susceptibility (panel (a), log scale), re¬ 
flection (panel (b), linear scale) and transmission (panel (c), 
linear scale) as functions of the reduced detuning of ®®Rb 
with A 0/7 = 21 . Vertical red dashed lines show the posi¬ 
tions of destructive interference between the different transi¬ 
tion dipoles. 


Overlapping of the transitions to F' = 2 and F' = 3 
from the ground states result in two Fano profiles in the 
susceptibility. Two sharp minima in the reflection spec¬ 
trum are the signatures of destructive interferences be¬ 
tween the transition dipoles. Concurrently we see two 
narrow peaks in the transmission spectrum. These fre¬ 
quencies at which DIET takes place are shown as red 
dashed lines in the Figure. 

Further increase of the dipole-dipole interaction mixes 
the F' = 2 and F' = 3 excited states into a single excita¬ 
tion state which is coupled to two ground states (F = 2 
and F = 3 of 5 ^S'i/ 2). It leads to DIET due to the over¬ 
lapping of transitions to the excited state from two differ¬ 


ent ground states. Eig. 9 shows such an effect in strong 
interacting samples of ®®Rb with A 0/7 = 200. Panel (a) 
in Fig. 9 is the modulus of the susceptibility, panel (b) 
and (c) are the reflection spectra for slabs of thicknesses 
i = 600 nm and £ = 2.5 /rm, respectively, and panels (d) 
and (e) are the corresponding transmissions. 



Reduced Detuning 5 


FIG. 9. (Color online) Susceptibility (panel (a), log scale) of 
a slab of ®®Rb atoms with A 0/7 = 200. Reflection spectra 
(panel (b) and c) for thickness I = 600 nm and I = 2.5 pm 
and corresponding transmission spectra in panel (d) and (e) 
as functions of the reduced detuning. A zoom of the suscep¬ 
tibility is given in the inset of panel (a) for better view of the 
faint Fano profiles associated with the holes in the reflection 
spectra. 


The susceptibility shows one dominant and two faint 
Fano type profiles. The latter clearly shown in the in¬ 
set in panel (a) are due to the splitting of the excited 
states. The dominant Fano profile corresponds to the 
overlapping of two coupled dipoles F = 3 —>■ F' = 2, 3 
and F = 2 —>■ F' = 2, 3. It leads to DIET which is well 
resolved in the reflection (Eig. 9(b) and (c)) and trans¬ 
mission (Fig. 9(d) and (e)) spectra of slabs of thickness 
600 nm and 2.5 pm. This can be understood if we assume 
that the DI transitions are a mixture of two three-level 
systems having the same excited states. In such a picture, 
the dominant DIET in the strong interaction regime can 
be interpreted as an interference effect between the tran¬ 
sitions to a single excited state from two different ground 
states. The two small Fano profiles in the susceptibility 
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are due to the splitting of the excited states. They also 
result in DIET which explains two small transmission 
peaks on the left and right of the dominant DIET signal 
from the 600 nm slab. These two transmission peaks are 
not resolved for the 2.5 /im slab due to damping effects. 
The system, therefore, acts almost as if it was comprised 
of three-level emitters [6]. However the splitting of the 
excited state leaves spectral holes in the reflection spectra 
(Fig. 9(c)). 

V. POTENTIAL APPLICATIONS 

The destructive interference of two overlapping tran¬ 
sitions can be used to modify the spectral envelope of 
an incident pulse, i.e. as a pulse shaper. If we consider 
a slab of interacting multi-level quantum emitters with 
a reflection window wider than the FWHM, 5oj of the 
incident pulse, the collective response of the quantum 
emitters induce a selective reflection and transmission of 
the incident pulse. This is illustrated in Fig. 10 for the 
D1 transitions of ®^Rb at Ap/q = 200, with a pulse of 
FWHM, 5u) = Att GHz. Panel (a) in Fig. 10 compares 
the normalized reflected pulse envelope (red line) with 
the incident pulse shape (blue dashed line) of such a pulse 
with the central wavelength 794.98 nm for a 600 nm thick 
slab of ®®Rb. Panel (b) is the same as panel (a) but for 
a 2.5 /rm thick slab. The normalized transmitted pulse 
envelopes (red line) are shown in panel (c) for 600 nm 
and in panel (d) for 2.5 /rm. A zoom around the central 
wavelength is given in the inset of panel (d) to see the 
shape of the pulse transmitted though the 2.5 p,m slab. 

From the Figs. 10(a) and (b) it is clear that the re¬ 
flection from the 600 nm slab is not as efficient as the 
reflection from the 2.5 /xm slab since this slab transmits 
a large part of the incident pulse (Fig. 10(c)). It is due 
to the mismatching between the slab thickness ^ = 600 
nm and the transition wavelength Aq — 795 nm that cre¬ 
ates long tails of Fabry-Perot modes near the reflection 
window which modify the transmission probability (see 
Fig. 9(b) and (d)). But for the 2.5 ^m slab, the con¬ 
dition £ > Ao is satisfied leading to the wide reflection 
window. This minimizes the transmission of the pulse of 
width 5uj < 3Ao through the slab except at the positions 
of DIET. 

If the system is comprised of emitters with many opti¬ 
cally active excited states that can be populated by the 
incident laser pulse (such as molecules with ro-vibrational 
levels for instance [31]), each transition interacts with 
the others and interferes constructively or destructively 
at different frequencies leading to reflected pulses with 
more than a single frequency amplified or removed from 
the spectrum. In particular, spectra of the the transmit¬ 
ted pulse in Fig. 10(c) and (d) show that DIET can vary 
dramatically depending on the system parameters. 


In systems comprised of interacting multilevel emitters 


I = 600 nm i = 2.5 /im 



Reduced Detuning (5 Reduced Detuning (5 


FIG. 10. (Color online)Reflected (panels (a) and (b)) and 
transmitted (panels (c) and (d)) pnlse shapes from the layer 
of ®®Rb atoms at Ao/q = 200 as functions of the reduced 
detuning &. The blue dashed lines in the panels are the nor¬ 
malized incident pulse and the red lines are the reflected and 
transmitted pulse envelopes. Panels (a) and (c) are for a slab 
of £ = 600 nm and panels (b) and (d) are for £ = 2.5 fj,m. 


we observe a minimum in the reflection window with a 
highly nonlinear dispersion (see Fig. 8 and 9). It opens 
several intriguing applications, including slow light [20, 
44]. Fig. 11 shows the slowing down of the transmitted 
pulse in a layer of ®^Rb with Ap/q = 200. Panel (a) 
shows the refractive index. The blue solid line is the 
imaginary part and the red dashed line is the real part 
of the refractive index. Panel (b) shows the group index 
ng{uj) = Re[n(a;) -I- uj{dn/du})] and the group velocity 
Vg(u}) = c/ng{u}) is shown in panel (c). 

The system has a very large group index Ug within the 
window of frequency where DIET occurs (see the Fig. 
11(d) and the inset of Fig. 10(b)). It peaks at 1.55 x 10^ 
similar to a recent experiment on slow light [44] and the 
group velocity (panel (c)) drops below 20 m/s. But in 
the system we consider this lowest velocity may not be 
observed since the imaginary part of the refractive index 
(see panel (a) in Fig. 11) is large that set the system as a 
lossy medium of light at this frequency. It also explains 
the quick drop in the transmission pulse shown in the in¬ 
set of Fig. 10(d). The transmitted pulse peaked around 
the reduced detuning (5 = 13 dies out at <5 = 20. Within 
this short window of the reduced detuning the group in¬ 
dex Ug is of the order of 10® and the velocity is decreased 
down to 40 m/s. Within this window, the imaginary part 
of the refractive index is small enough to allow the exper¬ 
imental observation of slowed down transmitted pulse. 
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Reduced Detuning 5 

FIG. 11. (Color online) Slow light effect associated with DIET 
in D1 transitions of ®®Rb for Aq /7 = 200. Panel (a) shows the 
refractive index: the solid blue line is the imaginary part and 
the dashed red line is the real part, panel (b) shows the group 
index of the layer and panel (c) shows the group velocity as 
functions of the reduced detuning S. 


VI. CONCLUSION 


When damping is sufficiently low, as in atomic or 
molecular vapors at low temperatures, for instance, nar¬ 
row transmission windows are observed. Such narrow 
windows in an otherwise completely opaque material are 
due to quantum interferences between different dipoles 
corresponding to different induced transitions. A clear 
evidence of such an interaction is the presence of Fano 
profiles in the susceptibility of the system. It is shown 
that this effect is amplified when the the dipole-dipole 
coupling increases. We call the transparency of the sys¬ 
tem induced by the coupling between different kinds 
of dipoles Dipole Induced Electromagnetic Transparency 
(DIET). We note that DIET is similar to the well-know 
electromagnetic induced transparency (EIT) [40, 41] the 
former however is fundamentally different from EIT. EIT 
requires a strong laser pulse that couples two different 
quantum states, while DIET is inherently internal, i.e. 
requires overlapping resonances induced by strong dipole- 
dipole couplings. 

The presence of additional levels such as ro-vibrational 
levels in molecules for instance or multiple optically ac¬ 
tive electronic states in case of atoms can change the 
response of the system due to the presence of additional 
energetically allowed transition dipoles which leave their 
signatures as spectral holes in the reflection spectrum 
and as narrow transmission peaks in the transmission 
spectrum. Such windows are induced via coherent can¬ 
cellation of the dipoles through a destructive quantum 
interference effect. In addition, we have shown that this 
effect can be controlled by changing the number density 
of quantum emitters. This destructive interference of two 
or more transition dipoles can be used for shaping laser 
pulses or for slow light. 


In this manuscript, we study collective effects in sys¬ 
tems of interacting multilevel quantum emitters and their 
dependence on different physical parameters such as den¬ 
sity. Such effects are studied both numerically using 
Maxwell-Bloch equations and analytically using a semi- 
classical model. 
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